Collisional oscillations of trapped boson-fermion mixtures approaching collapse 
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We study the collective modes of a confined gaseous cloud of bosons and fermions with mutual 
attractive interactions at zero temperature. The cloud consists of a Bose-Einstein condensate and a 
spin-polarized Fermi gas inside a spherical harmonic trap and the coupling between the two species is 
varied by increasing either the magnitude of the interspecies s-wave scattering length or the number 
of bosons. The mode frequencies are obtained in the collisional regime by solving the equations of 
generalized hydrodynamics and are compared with the spectra calculated in the collisionless regime 
within a random-phase approximation. We find that, as the mixture is driven towards the collapse 
instability, the frequencies of the modes of fermionic origin show a blue shift which can become very 
significant for large numbers of bosons. Instead the modes of bosonic origin show a softening, which 
becomes most pronounced in the very proximity of collapse. Explicit illustrations of these trends 
are given for the monopolar spectra, but similar trends are found for the dipolar and quadrupolar 
spectra except for the surface (n = 0) modes which are essentially unaffected by the interactions. 

PACS numbers: 03.75.Kk, 03.75.Ss, 67.60.-g 



I. INTRODUCTION 

Boson-fermion mixtures have been ^produced with 
atomic gases in several experiments 0, 0, H 0, la, E| and 
the deep quantum-degeneracy regime has been reached 
for both the bosonic and the fermionic component 0, @ • 
The aim of fermion cooling in these mixtures has been the 
realization of a superfluid state |3, |lfj . Even though the 
mixtures are very dilute, the interactions between bosons 
and fermions play important roles and in particular the 
s-wave collisions between the two species are exploited 
in the evaporative cooling process to obtain rethermal- 
ization. Mean-field interactions affect the expansion rate 
of the cloud and its thermodynamic properties 0, ^| . 

For a 87 Rb- 40 K mixture the boson-fermion interactions 
are large and attractive [12j , while in other systems they 
could be tuned to become attractive by exploiting Fesh- 
bach resonances Q, ^| . Boson-fermion attractions 
enhance the overlap between the two species, favoring in 
principle a boson- mediated fermion pairing [la . fl7| , and 
ultimately lead to the collapse of the cloud as the boson- 
fermion coupling strength overcomes the Fermi pressure. 
Collapse has been observed experimentally in the 87 Rb- 
°K mixture as a sudden disappearance of the fermion 
cloud when the number of bosons, and hence the cou- 
pling between the two species, is increased 01- The 
equilibrium properties and the phase diagram of a mix- 
ture with attractive interactions have been investigated 
by several authors 0, |2(], l2ll , finding good agree- 
ment between the experimental data and a mean-field 
theoretical description [23ll24l| . 

The study of the collective modes of a boson-fermion 
mixture with repulsive interactions has proved to be a 
good indicator of the quantum transition to a spatially 
demixed state |25j . For the case of attractive interactions 
some analyses using semi-analytical methods [2(| H?) 
have suggested that the lowest-lying monopolar mode 
should show a softening in the approach to collapse. A 



similar frequency softening is found for a pur e condensate 
with increasingly large boson attractions |28ll29| . In con- 
trast, the numerical solution of the equations of motion 
for boson-fermion mixtures in the collisionless Random- 
Phase- Approximation (RPA) does not give indications of 
frequency softening for the choice of parameters adopted 
in the calculations 30]. Moreover, it is found that when 
the number of bosons is much larger than that of fermions 
some families of modes are strongly blue-shifted as a con- 
sequence of the increase of the particle densities in the 
cloud approaching collapse |3l|. However, the RPA cal- 
culations could not be brought to the very proximity of 
collapse because of very severe numerical difficulties. 

In this paper we study the collective modes of a 
trapped boson-fermion mixture with attractive interac- 
tions using the equations of generalized hydrodynamics. 
These allow us to span the entire range of boson-fermion 
interaction strength up to the collapse point. Although 
the hydrodynamic equations are strictly valid only in the 
collisional regime, their predictions also serve as a guide- 
line to interpret the RPA spectra. By evaluating the 
frequencies of the first ten low-lying modes we find that 
some modes undergo a blue shift which can become quite 
large at large numbers of bosons, while some other modes 
display a softening which becomes pronounced in a very 
narrow window in the proximity to collapse. 



II. EQUATIONS OF GENERALIZED 
HYDRODYNAMICS 



We consider a model for a dilute mixture of two species 
of alkali atoms inside harmonic traps, one species being 
spin-polarized fermions of mass mj? and the other bosons 
of mass tub- The boson-boson and boson-fermion in- 
teractions are described by contact potentials involving 
the coupling constants Qbb = 47rfi, 2 clbb / rn-B and Qbf = 
2irh 2 asF /m r , where qbb and qbf are the corresponding 
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s-wave scattering lengths and m r = msmp / '(mp +ms). 
The fermion-fermion interactions are negligible, since the 
Pauli principle forbids collisions in the s-wave channel 
between fermions of the same spin. In the following we 
treat the case of boson- fermion attractions (clbf < 0) 
and boson-boson repulsions (clbb > 0). 

The dynamics of the mixture in the collisional regime is 
described by the equations of generalized hydrodynamics 
as in Ref. |25|. These are (i) the equations for particle 
conservation, 



dtPa + V • j CT = 0, 



(1) 



relating the particle density p a (r,t) to the current den- 
sity j CT (r,i) for each component (a = B,F); and (ii) the 
quantum Navier-Stokes equations for momentum conser- 
vation, which arc 



m a dtia = Pa (F a ~ W CT - g B F Vpa) 



(2) 



to linear order in the velocity fields. Here a denotes 
the component different from a, V^-(r) = m rT u 2 r 2 /2 
are isotropic trapping potentials, and F cr (r, t) are self- 
consistent internal forces. We take for the latter the ex- 
pressions 



F B = -V 



and 
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(4) 



with A = h 2 (67t 2 ) 2 / 3 /2tof- The equations of general- 
ized hydrodynamics are thereby closed by means of a 
local-density approximation on the kinetic stress ten- 
sors, which also includes surface kinetic contributions 
through the last term on the rhs of Eqs. © and 10} 25] . 
These surface correction terms transcend the standard 
Thomas-Fermi approximation and for bosons the present 
approach is in fact equivalent to the full time-dependent 
Gross-Pitaevskii equation. 

The equilibrium state of the mixture is given by the 
stationary solutions of Eq. J5J. These also correspond to 
a minimum of the energy functional 



E[pb,Pf] 



d 6 r 
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(5) 



where the surface kinetic-energy terms are £g = 
fr 2 \Vy/pB\ 2 /2m B and £ F = h 2 \ V^p^\ 2 /6m F . 

At increasingly large boson-fermion attractions the 
densities of the two species increase in their overlap re- 
gion and collapse occurs when these attractions overcome 



the Fermi kinetic pressure and the boson-boson repul- 
sions. Collapse is identified in the numerical search for 
the equilibrium state as the point where it is no longer 
possible to find a stable minimum for the energy func- 
tional in Eq. (JSJ). We have verified that the location of 
collapse found in this way is well approximated by the 
estimate Il9ll32l 



L BF I 



&BB 

a kp 
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(G) 



where kp = (6n 2 pi?(0)) 1 ^ 3 is the Fermi wave number 
estimated from the fermion density at the center of the 
trap and a = [3 1/3 /(27r) 2 / 3 ](r7i F + TO B ) 2 /(4m F TO B ). 

In the following we characterize the approach of the 
mixture to the collapse instability by following the be- 
havior of its low- lying collective oscillation frequencies. 
The normal modes of Eqs. £Q) and (0 are found by ex- 
panding yO CT (r,t) around the equilibrium state p^ir) as 
Pcr(r,t) = p"(r) + Sp a (r) e 1 " 1 , linearizing Eqs. (Q) and 
@, and Fourier transforming with respect to the time 
variable. This yields the coupled eigenvalue equations 

m a Lu 2 6p a = V • 0° SF a ) - 9 bfV- (p° V<5p 5 ), (7) 

where 8F a are the forces obtained as 8F a = F CT [p° + 
Sp a ] — F a [p°] to linear order in the density fluctuations 
(for their expressions see Ref. [2{|). We also determine 
the solutions of the eigenvalue equations in the case where 
the dynamical coupling between the two components is 
neglected. Hereafter we shall call these solutions the un- 
coupled boson and fermion modes and calculate them 
by solving Eq. Q) for a = B and F with the second 
term in the rhs set to zero. This calculation will help 
us in attributing bosonic or fermionic character to the 
coupled modes of the mixture by comparison with the 
frequencies of the uncoupled modes. Even though the 
bosonic and fermionic oscillations that are obtained in 
this way are dynamically uncoupled, the effect of the 
mutual boson-fermion interaction is included in their cal- 
culation through the use of the density profiles of the 
mixture at equilibrium. 

The numerical procedure that we have used can be 
summarized as follows: (i) we find the equilibrium den- 
sity profiles by minimizing the energy functional JSJ using 
the steepest-descent method; (ii) we project Eq. J7J) into 
subspaces of different angular momentum / by factorizing 
the density fluctuations as Sp a (r) = Sp l a (r) Yi m (f) where 
Yi m (f) are the spherical harmonics; and (in) we solve 
the coupled equations in a given Z-subspace by means of 
standard linear-algebra routines 33] . 



III. COLLECTIVE MODES 



We have solved Eq. (J7J) for two sets of experimentally 
relevant system parameters, corresponding to a 7 Li- 6 Li 



mixture |^| and to a 87 Rb- 40 K mixture [18J. In both 
cases we have set the frequency of the isotropic trap at 
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the geometric average of the frequencies in the experi- 
mental setup and have adopted the experimental value of 
the boson-boson scattering length. In our calculations we 
follow two different routes to reach collapse. We drive the 
7 Li- 6 Li mixture to collapse by varying the value of o,bf 
from positive to strongly negative, as could be achieved 
by means of a Feshbach resonance. In the 87 Rb- 40 K mix- 
ture instead we keep clbf fixed and increase the number 
Nb of bosons, as is done in the experiments of Modugno 
et al. 18]. 

A. The 7 Li- 6 Li mixture 

Figure ^ shows the frequencies of the low-lying 
monopole (I — 0) modes in the 7 Li- 6 Li mixture as func- 
tions of aBF, with the choice of parameters lvf = lob = 
2ir x 1000 s -1 , cibb = 0.27 nm, and particle numbers 
N F = 10 4 and N B = 10 6 . Pronounced blue shifts are ob- 
served at intermediate values of the boson-fermion scat- 
tering length both in a set of cigcnfrcquencies of Eq. Q 
(circles) and in the frequencies of the uncoupled fermionic 
modes (dashed lines). This effect is a consequence of 
the increase of fermion density in the central part of the 
trap as the strength of the coupling is increased. For 
our choice of parameters the bosonic density profile is al- 
most unaffected in this range of boson-fermion coupling 
and essentially acts as an effective attractive well for the 
fcrmions. The proximity of these modes to the uncou- 
pled fermionic modes indicates that the effect of the dy- 
namical coupling is negligible, thus confirming the inter- 
pretation in terms of static mean-field effects. At larger 
values of \cibf\ the other eigenfrequencies of Eq. 10 (cir- 
cles) show an increasing departure from the uncoupled 
bosonic modes (solid lines) and tend to rapidly soften as 
collapse is approached. We conclude that, in contrast 
to the blue-shift of the fermionic modes, the softening of 
"bosonic" modes is a truly dynamical signature of the 
impending collapse. 

Similar trends of the mode frequencies are also found 
for the dipole (I = 1) and quadrupole (I = 2) oscilla- 
tions, except that the lowest (n = 0) mode is in both 
cases essentially unaffected by the interactions. For the 
dipole mode this behavior is an exact consequence of the 
generalized Kohn theorem [34l [35J, while for modes of 
higher I it is strictly valid only in the Thomas-Fermi 
limit. This can be explicitly verified from Eq. by 
making the Ansatz 8p a (r) oc r l Yi m (f) for the density 
fluctuations and by exploiting the property ~S/ 2 5p a = 
for these modes. 



B. The 87 Rb- 40 K mixture 

We evaluate the collective modes of a 87 Rb- 40 K mix- 
ture with scattering lengths cibb — 5.5 nm and <xbf — 
— 21.7nm. Although both components are inside the 
same magnetic trap, their trapping frequencies differ 



considerably as a consequence of the large difference in 
atomic masses. In the numerical calculations we take 
lub = 27r x 90.9 s _1 and uj f = 27r x 134s -1 . Collapse is 
approached by varying Nb from 4 x 10 3 to approximately 
8 x 10 4 , while keeping N F = 2 x 10 4 fixed. 

In Fig. [21 we show the frequencies of the low-lying 
monopolar modes for this mixture. At increasing val- 
ues of the boson number we find that the frequencies of 
the bosonic modes, identified as those that are closer to 
the uncoupled bosonic fluctuations, show a monotonic 
decrease which becomes very pronounced very close to 
collapse. On the other hand, the fermionic modes only 
show a modest blue shift. This is not as large as in our 
results for the 7 Li- 6 Li mixture, since the numbers of par- 
ticles of the two species are here more similar. 

IV. SPECTRAL FUNCTIONS 

We have also calculated the whole spectral functions 
under the effect of the perturbing fields 5U a (r)e lult , in or- 
der to make a comparison of the collisional spectra with 
those of the mixtures in the RPA collisionless regime. 
Peaks in the spectral functions indicate the location of 
collective modes, while their height measures the prob- 
ability of a transition from the equilibrium state to an 
excited state with energy Huj. 

In the collisional regime the spectral peaks are located 
in correspondence of the eigenfrequencies loi of Eq. J7J, 
while the strength of a transition at a frequency to is esti- 
mated from the corresponding density fluctuations 5p l a (r) 
as 

S^}(u)=-hmJ2^f^ I d\5U;{v)8pl,{v) (8) 

i 1 

(see Appendix). Here the sum is performed over the 
whole spectrum of Eq. (Q and Z l a is given by 

Z F = I d 3 r (p F f 1/3 Sp F V ■ (p F WSU F ) (9) 

6m F J 

and 

Z l B = — I d 3 r Sp l B V • (p° B VSU B ) ■ (10) 
m B J 

Since we are using in Eq. © the equilibrium profiles 
and the mode frequencies of the fully coupled dynamical 
equations, we expect it to yield the spectral functions 
quite accurately. 

More generally, the spectral functions are related 
to the density-density responses x<r<r'( r > r '> w ) by tne 
fluctuation-dissipation theorem, 

S aa >(w) = -~Im J d 3 rd 3 r' 6U;(r) X aa>(r,r',Lu)5U a ,(r'). 

(11) 

We evaluate the response functions in the collisionless 
regime within the RPA, following the approach described 
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in Refs. 0,|3^|. For a monopolar drive (<5i7 CT (r) oc r 2 ) 
we calculate both the fermionic response S l fp(ui) and the 
bosonic response S l jfg(uj), corresponding to excitations 
induced by applying the drive separately on the fermions 
or on the bosons. 

In Figs. and 0] we report the spectral functions of 
the 7 Li- 6 Li mixture as obtained from the RPA formalism 
(top panels) and from Eq. JSJ for the collisional regime 
(bottom panels) at various values of cibf an d with the 
same choice of parameters as in Fig. ^ hi the absence of 
boson-fermion coupling a fermionic perturbing field ex- 
cites the pure n = mode of the fermions, which gives 
the only contribution to the FF spectra. In contrast, 
the BB spectra contain in the same limit several weak 
peaks at higher frequencies above the pure n = mode 
of the bosons, owing to the boson-boson interactions. For 
a BF 7^ the RPA spectra are very complex, from the lift- 
ing of degeneracy in the discrete-level structure which is 
due to the interactions with the bosons. The fermionic 
peaks can be recognized as broad fragmented contribu- 
tions of largest height in the fermionic spectral functions, 
while the bosonic peaks are usually isolated and are most 
marked in the bosonic spectra [3 1| . 

Quite interestingly, the spectra in the two regimes dis- 
play similar gross features. At small values of \clbf\ 
the peaks are located near the frequencies of the higher- 
order monopolar modes, but at larger values they tend 
to spread over the whole range of frequency. This phe- 
nomenon occurs for the same choice of parameters for 
which the blue shift of the collisional fermionic modes 
becomes pronounced and several mode crossings occur. 
For aBF — — 3 nm the first fermionic peak is strongly 
blue-shifted while the low-frequency part of the spectrum 
is dominated by bosonic modes having considerable os- 
cillator strength in both bosonic and fermionic spectral 
functions. This indicates a strong dynamical coupling of 
the density fluctuations of the two species in a regime of 
parameters where the uncoupled frequencies are close to 
the coupled solutions for the collisional fermionic modes 
(see again Fig. ^| . A careful analysis also shows that for 
aBF = — 3 nm the frequencies of the bosonic modes in 
the RPA spectra are slightly red-shifted, showing the first 
indications of a trend towards softening on the approach 
to collapse. 

Finally, in Fig. [5] we compare the collisional results 
for the monopolar spectral function (bottom panel) with 
those of the RPA (top panel) for a 87 Rb- 40 K mixture 
with the same choice of parameters as in Fig. [21 The 
RPA spectra show a fragmented fermionic contribution 
around the frequency 2ujp of the bare monopolar oscil- 
lation as well as two main bosonic peaks at its sides. 
The frequency of these two peaks tends to decrease with 
increasing similarly to what is found for the same 
modes in the collisional regime. This suggests that such 
frequency softening is a signature of the incipient col- 
lapse. 



V. SUMMARY AND CONCLUDING REMARKS 

In summary, we have used the equations of generalized 
hydrodynamics to study the collective modes of trapped 
boson-fermion mixtures with mutual attractive interac- 
tions driving the mixture towards the collapse instability. 
We have focused on two specific systems of experimental 
interest, for which we have chosen two different routes to 
collapse. 

In both cases we find that the collective spectra show 
a frequency softening of a family of modes of bosonic 
nature as a signature of the incipient collapse. This soft- 
ening becomes most pronounced in a very narrow region 
of parameters near collapse. A second effect of increas- 
ing the attractive coupling between the two species is a 
blue shift of the modes of fermionic origin, which becomes 
more evident for large numbers of bosons and reflects the 
compression of the fermionic cloud from the interactions 
with the bosons. 

A comparison of the hydrodynamic spectra with the 
spectra calculated in the collisionless regime within the 
random-phase approximation suggests that both the blue 
shift of the "fermionic" modes and the softening of the 
"bosonic" modes are general features of the dynamics of 
the mixture. By comparing the height of the peaks in the 
bosonic and fermionic spectral functions we also conclude 
that the two families of modes can be most efficiently 
excited by applying the driving fields separately on the 
species. 

Our analysis has been restricted to the linear regime, 
neglecting nonlinear and beyond-mean-field effects which 
may start to play a role close to collapse. We hope to 
address these issues in future work. 
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APPENDIX: SPECTRAL FUNCTIONS IN THE 
COLLISIONAL REGIME 

The spectral functions of a mixture in the collisional 
regime can be obtained by adding the perturbing fields 
SU a (r)e luJt to the trapping potential V a (r) in Eqs. JTJ and 
P|l and then linearizing the equations of motion in terms 
of the induced density fluctuation 5p a (r,uj). This leads 
in Fourier transform with respect to time to the equation 

m a uj 2 5 Pa = V-(p° SF a )-g BF V-(p° V<5p ff )-V-(p° V<ft/ CT ) 

(A.I) 

for each component of the mixture. Here, at variance 
from Eq. J7J uj is fixed by the frequency of the drive. 
Equation (jA.lfl can be solved by expanding (5p CT (r, w) in 
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the basis of the eigenvectors 8p % a (r) of Eq. J7J as 

6 Pa (r,u;) = 6pl(r). (A.2) 



In turn, the expansion coefhcients C 1 are obtained by 
projecting Eq. IjA.ljl into the subspace generated by a 
given fluctuation £p£.(r). This procedure requires defin- 
ing a scalar product between two eigenmodes of Eq. J7J , 
which we shall denote by the braket (5p l \5p'). The pro- 
jection yields 

C » . (VIV-W^)> (A .3) 



where uji is the frequency corresponding to the eigenvec- 
tor 8p % a (r). 

A generic scalar product can be written as 
(SpW) = E/ d 3 r6pl(v)w aa ,(r)6pi,(r) (A.4) 

ct, a' 

where w aa i (r) are suitable weights to be determined. 



Here we take for the scalar product (|A.4(1 that of a 
bosonic and of a fermionic cloud in the Thomas-Fermi 
limit with vanishing mutual interactions. We thus set 
Wff — 2A(p° F )^ 1 / 3 /3niF and wbb = 9BBl°n^B and ne- 
glect the off-diagonal weight functions. This choice yields 



CM 



ujf - UJ 2 



(A.5) 



for fermions and bosons, where Z % a is given in Eqs. © 
and (|10|l . We have taken 8p % G as normalized to unity ( 
(SpW) = 1). 

The spectral functions are then calculated from the 
definition 

£r<r'M = ~lm J d 3 rSUZ{r)6 Pl7 ,{r,u;) (A.6) 

where 8p a i is evaluated by setting SUgi = 0. The use of 
Eqs. (| A. 2|> and HA.5|) in Eq. I|A.6|) directly leads to Eq. 
(151) in the main text. 
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FIG. 1: Monopolar frequencies w (in units of the bare boson- 
trap frequency ujb) of the collisional collective modes as func- 
tions of the mutual scattering length clbf (in nm) for a 6 Li- 
7 Li mixture with N F = 10 4 and N B = 10 6 . The solid and 
dashed lines show the frequencies of the dynamically uncou- 
pled modes for bosons and fermions, respectively. The tri- 
angles correspond to the values of o,bf for which we have 
performed the RPA calculations of Figs. |3] and 0] 
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FIG. 2: Monopolar frequencies u) (in units of ujb) of the 
collisional collective modes as functions of the number Nb of 
bosons (in log scale) for the 87 Rb- 40 K mixture with Nj? = 
2 x 10 4 . The solid and dashed lines show the frequencies of 
the dynamically uncoupled modes for bosons and fermions, 
respectively. The triangles correspond to the values of Nb for 
which we have performed the RPA calculations of Fig. |S] 
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FIG. 3: Monopolar spectrum (in log scale and arbitrary 
units) for fermions as functions of uj (in units of ujb) in a 6 Li- 
7 Li mixture with Nf = 10 4 , Nb = 10 6 and for various values 
of clbf- The spectra are plotted for the sake of clarity with 
a spectral width of order 1CP 5 ujb- The vertical dotted lines 
indicate the values of the bosonic and fermionic monopolar 
modes in the absence of boson-fermion coupling. Top panel: 
spectral response within the RPA; the arrows indicate the 
frequencies of the collisional modes with non-negligible spec- 
tral weight for the same parameters. Bottom panel: spectral 
response in the collisional regime; the arrows denote all col- 
lisional modes for the same parameters with label B or F to 
indicate the type of mode as identified from its proximity to 
an uncoupled collisional mode (the labels BF indicate am- 
biguous cases). 
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FIG. 4: The same as in Fig.EJfor the bosonic response. 
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FIG. 5: Monopolar spectral response (in log scale and arbi- 
trary units) from the RPA for fermions (top) and bosons (bot- 
tom) as functions of w (in units of ujf or ujb) for a 87 Rb- 40 K 
mixture with clbb = 5.5 nm, clbf = — 21.7nm, Nf = 2 x 10 4 , 
and increasing values of Nb- The arrows indicate the colli- 
sional modes for the same parameters. The vertical dotted 
lines are the values of the bosonic and fermionic monopolar 
modes in the absence of boson-fermion coupling. 



